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Abstract. We present a group-theoretic criterion under which one may verify the Artin 
conjecture for some (non-monomial) Galois representations, up to finite height in the 
complex plane. In particular, the criterion applies to ^5 and representations. Under 
more general conditions, the technique allows for the possibility of verifying the Riemann 
hypothesis for Dedekind zeta functions of non-abelian extensions of Q. 

In addition, we discuss two methods for locating zeros of arbitrary L-functions. The 
first uses the explicit formula and techniques developed in BSOS^ for computing with 
trace formulae. The second method generalizes that of Turing for verifying the Riemann 
hypothesis. In order to apply it we develop a rigorous algorithm for computing general 
L-functions on the critical line via the Fast Fourier Transform. 

Finally, we present some numerical results testing Artin's conjecture for 6*5 representa- 
tions, and the Riemann hypothesis for Dedekind zeta functions of ^5 and A5 fields. 



1. Introduction 

1.1. Artin's conjecture. Let K/Q be a Galois extension and p : Gal(i^r/Q) — >• GL„(C) 
a non-trivial, irreducible representation of its Galois group. In |Art30j . Artin associated 
to this data an L-function L{s,p), defined initially for 3fJ(s) > 1, which he conjectured 
to continue to an entire function and satisfy a functional equation. By a theorem of 
Brauer |Bra47j . one now knows the meromorphic continuation and functional equation of 
Artin's L-functions. The question remains whether they can have poles in the critical strip 
< < 1. 

Artin established his conjecture for the monomial representations, those induced from 
a 1-dimensional representation of a subgroup; this of course includes all 1-dimensional p, 
in which case L{s,p) = L{s,x) for a Dirichlet character x- Although the conjecture has 
not been decided in any dimension > 2, more evidence is provided in dimension 2 by the 
Langlands-Tunnell theorem |Lan80tlTun81j . which affirms the conjecture for those represen- 
tations whose image in PGL2(C) is isomorphic to (tetrahedral) or 5*4 (octahedral); only 
the (icosahedral) case remains. When p is an odd icosahedral representation, meaning 
det p determines an odd Dirichlet character, infinitely many examples of Artin's conjecture 
are known from the work of Taylor et al. [BDSB ToH |Tay03| . 

Moreover, in the odd 2-dimensional case, there is an algorithm for verifying the con- 
jecture, as follows. By a construction of Deligne and Serre |DS74j . given a holomorphic 
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modular form / of weight 1, one may associate an odd 2-dimensional representation p such 
that L{s, f) = L{s,p). Conversely, every odd 2-dimensional p such that L{s,p) is entire 
arises from the Deligne-Serre construction. For any particular p, one can search for the 
associated form; once found, comparing the representation constructed by Deligne-Serre 
to p via an effective version of the Cebotarev density theorem allows one to deduce the 
conjecture for p. This and other related techniques have been carried out in a number of 
cases; see |Buh781 iKiSoil iJMnnl 16^1)2] . 

On the other hand, if one considers even 2-dimensional representations, the situation 
is somewhat different. There as well the conjecture has been established for all but the 
icosahedral cases. However, the correspondence is not with holomorphic forms, but rather 
Maass forms of eigenvalue | . Unfortunately, no analogue of the result of Deligne and Serre 
is known in that setting. Moreover, computation of the associated forms remains elusive; 
existing techniques (see e.g. |BSV05j ) only allow one to calculate Maass forms to within 
a prescribed precision, never exactly. Thus, at present this approach does not yield an 
algorithm for verifying Artin's conjecture. 

The apparent difference between these two cases leads naturally to the following question: 
Given a Galois representation p, is there an algorithm that will decide in finite time, with 
proof, whether L{s,p) is entire? Note that like the Riemann hypothesis, Artin's conjecture 
is falsifiable, i.e. it may be disproven by observing a counterexample, in this case a pole. 
The challenge is thus to find a way of demonstrating the conjecture when true. 

Although we are unable to provide a definitive answer to this question, one approach, at 
least for 2-dimensional representations, is suggested by a theorem from |Boo03bj : If a given 
2-dimensional p is not associated to a holomorphic or Maass form as above, then L{s,p) 
has infinitely many poles. In particular, once L{s,p) has at least one pole, it must have 
infinitely many. Unfortunately, the result is ineffective, in the sense that it does not predict 
where the first pole must occur. A natural question, therefore, is whether an effective 
version of this theorem exists. First, however, we must consider exactly what that would 
mean; since the only handle that we have on an Artin L-function in the critical strip is as 
the ratio of entire functions given by Brauer's theorem, it is not immediately clear that we 
can check its holomorphy at a zero of the denominator without a priori knowing a lower 
bound on the residue of any poles. 

In this paper we address precisely this issue, in Section |21 There we present a criterion 
which, when satisfied, yields an algorithm for verifying the holomorphy of an Artin L- 
function up to a given height in the critical strip. In particular, we give the first direct 
evidence (as far as we are aware) of holomorphy in the critical strip of an L-function 
for which the conjecture cannot be established through the methods mentioned above. 
Although our criterion is not always satisfied, we are in general able to deduce partial 
information, such as a bound on the multiplicities and residues of possible poles. Moreover, 
the limitations of the information that we obtain give an idea of the hypotheses that one 
would have to impose in any effective version of the converse theorem in order to make the 
above approach work. 
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1.2. Turing's method and the Riemann hypothesis. One application of our criterion 
is to the Riemann hypothesis for Dedekind zeta functions. Turing |Tur53j devised a method 
for checking the hypothesis in a bounded region for the Riemann ( function.^ The method 
depends on the simphcity of the zeros of (. Because of that, it is only directly extendable to 
Dedekind zeta functions of non-normal extensions of small degree (see |Tol97j ) or abelian 
extensions, for which it is more natural to verify the hypothesis for the associated Dirichlet 
L-functions instead (see jRum93j). 

Similarly, for a non-abelian extension one can factorize the zeta function into Artin L- 
functions of irreducible representations. As these are also expected to have simple zeros, 
Turing's method applies, provided one assumes the Artin conjecture. However, combining 
our criterion with Turing's method, we will in some cases be able to deduce the Riemann 
hypothesis and holomorphy of the relevant Artin L-functions simultaneously. In fact, as we 
will see, there are even cases where we may check the Riemann hypothesis without being 
able to verify Artin's conjecture. We carry out the necessary generalization of Turing's 
method in Section |3] 

1.3. Rigorous zero computations. In order to implement these ideas, we develop, in 
Sections El and El two methods of locating zeros of L-functions. The first uses the explicit 
formula and techniques developed for the Selberg trace formula in |BS05j . If one assumes 
the Riemann hypothesis, this method may be used with our criterion, in place of Turing's 
method, for verifying the Artin conjecture. More importantly, the explicit formula is clean 
to implement and yields estimates for low zeros quickly. It can thus serve as a check for 
later computations, or to fine tune the parameters of Turing's method for greater speed. 

The second method is a technique for fast, rigorous computations of L-functions on 
the critical line. This is a hard problem in general, basically because of the difficulty of 
providing uniform, effective bounds for the relevant Mellin transforms. By making use of 
the Fast Fourier Transform, our technique allows one to compute many values of the same 
L-function simultaneously, which is particularly appropriate for Turing's method. In doing 
so, we need only consider a single Mellin transform, making rigorous computation more 
practical. In addition, the method has complexity comparable to that of computing a single 
value by the approximate functional equation. 

Although our primary interest is in Artin L-functions, we carry out the details of Sections 
EllHandElfor arbitrary L-functions L(s), in the hope that the results may be useful outside 
of the present context. More precisely, we make the following assumptions, notations and 
conventions throughout: 

• L(s) is given by an Euler product of degree r: 



(1.1) 



Hs)= n 



1 



p prime 



(1 - Op^ip ^) ■ ■ ■ (1 - ap^rP 
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where the are complex parameters satisfying the individual bound \cepj\ < p^^'^, 
and the product is absolutely convergent for 3fJ(s) > 1. Further, for all but finitely 
many p, there is a pairing a a' such that lapjOp ,, | = 1. For the exceptional p, 
such a pairing exists for a subset of the apj, and those not in the subset satisfy 



Define 

r 

T^{s) := n-^/'T (0 , jis) := eN^^^"^^ J] ^^(^ + /^.), A(.) := 7(^)^(^), 

i=i 

where |e| = 1, is a positive integer and ^{fij) > — |. For a certain choice of these 
parameters, A(s) has meromorphic continuation to C, is a ratio of entire functions 
of order 1, and satisfies the functional equation 

A(s) = A(l-s), 



where for a complex function / we denote by f{s) the function f{s). Note that e here 
is the square root of the usual root number, and is only defined up to multiplication 
by ±1; we choose the value with argument in (— f , f ] ■ Including e as part of the 7 
factor makes A(s) real for 3fJ(s) = |, as can be seen from (jl.3j) . 
Let Q{s) be the analytic conductor. 



271 

Note that 7(5) satisfies the recurrence 7(5 + 2) = Q{s)'y{s). Further, we define 

7(1 - s) 



Xis) 



7(s) 



so that L{s) = — s). 

L[s) may have at most finitely many poles, which we assume to lie along the line 
3fJ(s) = 1. We label them 1 + with G zM, k = 1, . . . ,m, repeating with the 
appropriate multiplicity. Further, from the functional equation (jl.3p . each A^ will 
equal —fij for some j, counting multiplicity; in particular, m <r. We set 

m 

P{s) := J](s-A,), 

k=l 

SO that P{s)P{s — l)A(s) is entire. 

Some progress is known toward the Ramanujan conjecture for L; that is, there exists 
9 < ^ such that 

|«pj| < and > —9 

for all p, j. This assumption is not strictly necessary, as we could instead use average 
bounds of Rankin-Selberg type. However, bounds of the form ()1.7j) are now known 
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in the cases of greatest interest (automorphic L-functions |LRS99j ). and the results 
are easier to state and use assuming it. 

1.4. Numerical results. Finally, in SectionlHlwe describe the implementation of the above 
ideas and give some numerical results of tests of the Riemann hypothesis for a few and 
extensions in the region < 100. For the cases, this includes a verification of 

Artin's conjecture in the same region for the L-functions of all representations of the group. 



The idea for this paper arose from a conversation with Laurent Clozel, and it is a pleasure 
to acknowledge him here. I thank Harold Stark, who has considered problems of this nature 
in the past and whose interest motivated me to finish the project. Thanks also to Bob 
Griess, Martin Isaacs, Peter Sarnak and Kannan Soundararajan for helpful discussions. 



Let p : Gal(i^/Q) GL„(C) be a Galois representation, as in the introduction. Brauer's 
theorem expresses the L-function L{s,p) as a ratio N{s)/D{s), where N{s) and D{s) are 
Artin L-functions associated to sums of monomial representations. If f{s) is any holomor- 
phic Artin L-function, we have a formula for the number Nf{ti,t2) of zeros of / between 
heights ti and t2, from the argument principle: 



where C is the rectangular contour with vertices at 2 + iti, 2 + it2, —l+it2, —l+iti and 
counter-clockwise orientation. We also have available in this case algorithms to compute 
/ and /' at an arbitrary point in the complex plane; see Section |3| Thus, in principle we 
could compute (j2.H) exactly by numerical integration. Although (j2.H) has the advantage 
of applying in great generality, to do so would be inefficient and difficult to implement 
rigorously. In the special case that the zeros of / are simple, a much more efficient algorithm 
was given by Turing; see Section HI 

No matter how we arrive at the numbers Nf{ti,t2), there is always some uncertainty 
in the locations of the zeros of /. In fl2.1|) this is due to the fact that as ti approaches 
the ordinate of a zero, higher and higher precision is needed in order to compute /'// 
accurately. This is in line with the expectation that the general zero is transcendental, 
meaning that one can never know it exactly. 

For L(s, p), we can recover the net number of zeros (i.e. zeros minus poles) between 
heights ti and t2 as N]y{ti, ^2) — N£)(ti, ^2)- If Artin's conjecture is true then for every zero 
of D{s) there is a zero of N{s) at the same point. However, because of the uncertainty in 
the locations of the zeros of N{s) and D{s), from this computation alone we cannot rule 
out the possibility that L(s, p) has a pole with a zero very close by in the neighborhood of 
a zero of D{s). In other words, we can only observe the counts of net zeros in these small 
neighborhoods. 
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(2.1) 
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Fortunately, there are some restrictions on potential poles. For instance, the Dedekind 
zeta function of the extension, Ck{s), factors into Artin L-functions: 

(2.2) Cx(s) = n^(^'/^)''"'' 

p 

where the product is over all irreducible representations of Gal(-ft'/Q). Since Ck{s) is 
holomorphic (except for a simple pole at s = 1), we see that any pole of L{s,p) in the 
critical strip must be located at the zero of another function. More generally, if a is any 
representation, we have 

(2.3) L{s,a) = l[L{s,py^'''\ 

p 

where (■, ■) is the inner product on the space of characters, and by abuse of notation we 
write (cr, p) for (Trcr, Trp). When a is monomial, we again have L(s,cr) holomorphic with 
the possible exception of a pole at s = 1. 

This information is described most concisely by use of the Heilhronn (virtual) character. 
For So e C\ {1}, define 

(2.4) = oTds=soL{s, p) ■ Tr p, 

p 

where oids=SQL{s, p) := ReSs=s„^(s, p). Thus, 

(2.5) OTds=soL{s,(j) = {6sg,a) > for all monomial a. 

The study of Heilbronn characters leads to many useful results. For example, in |FM89] it 
is shown that 

(2.6) Y{oTds=soL{s,p)Y < [oTds=soCK{s)Y . 

p 

In particular, the zeros and poles of each L{s,p) are among the zeros of Ck{s). 

The idea now is to combine ()2.5|) with observations of net zeros. If we look in a small 
enough neighborhood of a zero of Ck{s), we expect to find one net zero for a single L{s,p) 
and no net zeros for the others. This is based on the assumption that the zeros of different 
irreducible Artin L-functions are distinct and simple. While such a statement is likely 
impossible to prove, we may use it as a working hypothesis to be tested at run time. This 
is analogous to assuming the simplicity of the zeros of ( in order to check the Riemann 
hypothesis. (Note that if there were a multiple zero of (, it is doubtful that one could 
distinguish it from a counterexample.) 

In other words, if the working hypothesis is true, then our net zero observations cor- 
respond to the character Tr p for some p. Thus, we have Tr p = ^<j^ + . . . + 6s„ , where 
Si, . . . ,Sn are the distinct zeros of Ck{s) in the neighborhood that we examine. We would 
like to conclude that there is just one such point, meaning that the actual zero counts agree 
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with our observations. Since the Heilbronn characters satisfy ()2.5|1 . it is enough to show 
that 

(2.7) Tr p 7^ Xi + X2 for virtual characters Xi 7^ with [xii c") > for all monomial a. 

The one notable exception to this philosophy is at the central point |, where there can 
be forced vanishing if p is self-dual (an example of which is given in |IArm72j ) . In that case, 
we expect one zero for each self-dual p with an odd functional equation, and no zeros for 
the rest. However, we can only determine the parity of the order of vanishing at |. This 
leads to the following replacement for condition ()2.7p at |: 

(2.8) Tr p 7^ Xi + 2x2 with Xi 7^ and (xi, o") > for all monomial a. 

p self— dual 
Hl-s,p)=-\{s,p) 

When ()2.7|) is satisfied for all irreducible representations p, we may check the holomorphy 
of all L{s,p) at any point at which the working hypothesis turns out to be true. We give 
a name to describe this situation: 

Definition 2.1. A finite group G is almost monomial ii, for each irreducible representation 
p, if Trp = Xi + X2 for virtual characters Xi such that (xi, o") > for all monomial cr, then 
either xi = or X2 = 0. 

The terminology is explained with aid of Figure 12.11 The plane represents the lattice of 
virtual characters, with the first quadrant being the monoid of characters, and the shaded 
cone the monoid generated by the monomial characters. We consider all virtual characters 
within 90 degrees of the cone, which in the figure is everything within the dashed lines. 
The group is almost monomial if this set is not much larger than the character monoid, in 
the precise sense that the irreducible representations, which are the coordinate axis vectors 
represented by thick arrows, remain indecomposable in this set. Equivalently, the monoid 
generated by the monomial characters should be close to the full character monoid. From 
the picture it is easy to see that any monomial group is almost monomial. 

One could argue that we should include condition ()2.8|) in our definition as well. We prefer 
to keep it separate, taking the view that it is more important to be able to demonstrate 
holomorphy at a generic zero of the denominator. Indeed, we have already seen that the 
L-function of a 2-dimensional representation cannot have a finite number of poles, so we 
do not lose much generality by excluding a single point. It is plausible that such a result 
holds for higher dimensions as well. Moreover, condition (j2.8|) seems usually to be weaker 
than almost monomiality; cf. Proposition 12.31 below. 

A potentially more serious issue is that N{s) and D{s) may have high order zeros at ^, 
in which case Turing's method does not apply. This could be remedied by computing the 
contour integral ()2.1|) around |, but we would like to avoid doing so. Fortunately, if the 
order of vanishing at | is at most 3, we can still conclude that we have the correct count 
by sign changes alone; that is because for a self-dual representation, if we miss a zero away 
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Figure 2.1. Lattice of virtual characters. The shaded cone is the monoid 
generated by monomial characters. 

from I then we must miss at least four such zeros. Fortunately again, in all cases that we 
consider, N{s) has at most three irreducible factors with a potential zero at |. 

Like monomiality, the notion of almost monomiality behaves well under some common 
group operations. In particular, we have the following. 

Proposition 2.2. If G is almost monomial then so are quotients of G and products G x H 
for any monomial group H. 

Proof. 1. Let be a normal subgroup of G and tt an irreducible representation of GjK. 
Suppose that Tryf = Xi + X2, with (xj, o-) > for all monomial u. Let tt, Xi be the lifts of 
71", Xi to G obtained by composition with the natural projection. Then tt is irreducible and 
Tr TT = Xi + X2- Further, if p is an irreducible representation of G then (xj, p) = unless p 
factors through GjK. If that is the case, let p denote the induced map on GjK. 

Now, \i a — Ind^A is a monomial representation then, by Frobenius reciprocity, we have 
(cr, p) — (Res^p, A) for all p factoring through GjK. If at least one of these is non-zero, i.e. 
A occurs in Res^p, then since p factors through G/K, \ must factor through H/H D K = 
HK/K. Let A denote the induced map on HK/K. Then (ResJ^^^p, A) = (Resgp,A>. 

Thus, a = Ind^'^^A satisfies (cr, p) = (a-,p). Therefore, {xi,(^) — (Xt^^) ^ 0- The 
conclusion follows by almost monomiality of G. 

2. Let pg and pn be irreducible representations of G and H, respectively, and suppose 
that Tt pg®Ph = (Tr pG)(Tr pn) = X1+X2 with {xi, cf)gxh > for all monomial cr. Taking 
the inner product over H with pn, we get Trpc = (xi, Ph)h + (X2, Ph)h- 
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Next, if ac is any monomial representation of G, we have {{Xi, Ph)h,<^g)q = {Xii^G ® 
Ph)gxh > 0, since ac ® Ph is monomial. Thus, since G is almost monomial, we have 
{Xii Ph) = for some i. 

Similarly, if p'^j^ is any other irreducible representation of H, we find = {xi^p'h)h + 
{X2,p'h)h- Thus, {xi,p'h) = {X2, p'h) = 0- Therefore Xi = for some i. □ 

The next proposition shows that the class of almost monomial groups is strictly larger 
than that of monomial groups. 

Proposition 2.3. The groups SL2(F3), and are almost monomial and satisfy ()2.8|) . 

Proof. These are shown with the aid of the computer algebra system GAP |GAP05j . We 
illustrate the general procedure for checking almost monomiality for a given group with 
the example ^5. Note first that has five irreducible representations, of dimensions 1, 
3, 3, 4 and 5. We use GAP to determine all monomial representations. In this case they 
are non-negative linear combinations of the vectors (1, 0, 0, 0, 0), (0, 0, 0, 0, 1), (1, 0, 0, 1, 0), 
(0,1,1,0,0), (0,1,0,1,1), (0,0,1,1,1), and (0,1,1,1,0), where the components indicate 
the multiplicities of the irreducible representations. We label the monomial representa- 
tions associated to these vectors ai, . . . , (J7. The first five form a Z-basis for the lattice of 
virtual characters, i.e. any virtual character x Kiay be written uniquely as an integral linear 
combination x = Yl^=i ^^iTr (Xj. 

Now, almost monomiality is equivalent to the assertion that for each irreducible repre- 
sentation p, whenever < {x, u) < {p, a) for all monomial a, we have either % = or 
X = Tr p. Using our integral basis, we investigate the solutions to 

5 

(2.9) 0< ^a;,(a„a,) < (p,a,) 

1=1 

for j = 1, . . . ,7. Restricting to j = 1, . . . , 5, we get an invertible system, i.e. the matrix A = 
[{ai, aj)) lies in SL^i^Z). We consider the vectors x = A'^y for all y = (yi, ... ,7/5) 
satisfying < yj < (p, dj). By construction, these satisfy (j2.9p for j = 1, . . . , 5. We check 
that the only x satisfying ()2.9p for j = 6, 7 are and A^^[{p, aj)) , corresponding to x = 
and X = Tr p, respectively. 

Similarly, for ()2.8|) we try all possible combinations of p having odd functional equation. 
We may exclude those whose L-functions may be expressed in terms of Dedekind zeta 
functions, for which the root number is always 1. For A^, the only non-trivial possibility is 
that the two 3-dimensional representations have odd functional equation. □ 

With the evidence provided by Propositions 12.21 and 12. 3[ one might hope that all groups 
are almost monomial. That is not the CcLSG, 8jS the counterexamples GL2(F3) and SL2(F5) 
show. SL2(F5) has irreducible representations of dimensions 1, 2, 2, 3, 3, 4, 4, 5 and 6, 
and it is the smallest group supporting an icosahedral representation (since A5 has no 2- 
dimensional representations), meaning that our criterion unfortunately does not apply to 
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checking the icosahedral case. In fact, one knows Artin's conjecture for all induced repre- 
sentations of this group; while they are not all monomial, the only exceptions come from 
a pair of tetrahedral representations, for which we have the Langlands-Tunnell theorem. 
Even with this added information, we cannot rule out the possibility of a simple pole with 
undetectably small residue at a zero of the L-function of the 6-dimensional representation. 
More precisely, we find with GAP that the induced representations are spanned by the 
twelve vectors 



The first seven of these are the monomial representations lifted from SL2(F5)/{±/} = A^, 
while the others give "new" information. One easily checks that for p the 6-dimensional rep- 
resentation, (j2.7p fails with Xi corresponding to any of the vectors (0, —1, 0, 0, 0, 0, 0, 0, 1), 
(0,0,-1,0,0,0,0,0,1) and (0,0,0,0,0,0,-1,0,1), i.e. the representations of dimension 2 
and one of dimension 4 can hide a pole at a zero of L{s,p). This shows in a strong sense 
that information from induced representations is in general insufficient to show Artin's 
conjecture. 

However, all is not lost concerning icosahedral representations. For a given icosahedral 
p, the adjoint square Ad(p) is a 3-dimensional representation with image isomorphic to 
A^. A result of Flicker |Fli94j implies that modularity of p is equivalent to that of Ad(p). 
(In fact, modularity of all representations of the underlying group follows from that of 
Ad(p) and its Galois conjugate, by known cases of functoriality; see |Wan03j .) Combining 
this fact with the GL(3) converse theorem, one could give a converse theorem for GL(2) 
using analytic properties of L{s, Ad(p) ® x) fo^' Dirichlet characters x- Weissman, in his 
undergraduate thesis |Wei99j , used this idea to give indirect evidence for the modularity of 
an even icosahedral representation. By Propositions 12 . 21 and 12 . 31 we see that in principle we 
may directly verify the holomorphy of these L-functions up to finite height. Moreover, an 
"effective" version of the GL(3) converse theorem (requiring, say, meromorphy of all twists 
and holomorphy of a finite number in a bounded region) would suffice to give an algorithm 
for verifying the conjecture in the icosahedral case. 

Unfortunately, there is the more practical problem that totally real A^ fields (those that 
give rise to even icosahedral representations) are very rare; the smallest known discriminant 
is far too large to test with current computers. Thus, for the A^ examples that we consider 
in Section El the Artin conjecture is already known. To test our criterion, we consider 
instead some examples of 5*5 extensions, which exist in much greater abundance. 

Finally, we note that in the course of verifying Artin's conjecture, the information that 
we collect implies that the zeros of each L{s,p) are simple and lie on the line 3fJ(s) = |. 
Thus, in the process we also verify the Riemann hypothesis for Ck{s)- Interestingly, we 
do not need to establish the holomorphy of all L{s,p) in order to do this; it is enough. 



(2.10) 



(1,0,0,0,0,0,0,0,0) 
(0,0,0,1,1,0,0,0,0) 
(0,0,0,1,1,1,0,0,0) 
(0,1,0,0,0,0,1,0,1) 



(0,0,0,0,0,0,0,1,0) 
(0,0,0,1,0,1,0,1,0) 
(0,0,0,0,0,0,0,0,1) 
(0,0,1,0,0,0,1,0,1) 



(1,0,0,0,0,1,0,0,0) 
(0,0,0,0,1,1,0,1,0) 
(0,0,0,0,0,0,1,0,1) 
(0,1,1,0,0,0,0,0,1). 
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for example, that they have at most simple poles. More precisely, in order to check the 
Riemann hypothesis around a generic zero of Ck{s) we need to have 

(2.11) Tr p 7^ Xi + 2x2 for virtual characters 7^ with (xi, o") > for all monomial a, 

which is a weaker condition than almost monomiality. In particular, we may still check the 
Riemann hypothesis for SL2(F5) extensions. 



3. Locating zeros via the explicit formula 

Let notation be as in the introduction, and define numbers c„ by —j;{s) = c^n"**, 
i.e. Cn = (logp) X]j=i '^pj n = a. prime power, and c„ = otherwise. Further, 
we enumerate the zeros of A(s) as p„ = | + 27„ for n G 7j, repeated with multiplicity. 
Weil's explicit formula relates the sequences {c„} and {jn}- Precisely, let g G C^(M) 
be a different iable function of compact support such that its Fourier transform h{t) : = 
I^oo di^)^^^^ dx is real for t G M. Then 



(3.1) 



neZ k=l ^ ^ 



2 + 



9(0)logiV + 2S 



n] 



This follows from the Cauchy integral formula and the functional equation; see 
Note that all terms of the formula may be put in terms of g; in particular, 

1 poo p/ /I \ 1 poo 

(3.2) -j ^[2+^' + 't)h{t)dt = -j^ \og{7Te^e'^-l))d{g{x)e-^'/'^^^^). 

This form is convenient for computation, since g has compact support. 

The important thing to note is that given a list of the c„ for n < , the explicit 
formula gives a method for evaluating ^(7n) for essentially any function h whose Fourier 
transform is supported in [— X, X]. When X is large, we may choose h to be narrowly 
concentrated around any particular point, and thus resolve features of the spectrum in 
places where the density of zeros is not too large compared to X; a variant of this technique, 
with explicit test functions (not of compact support), was worked out by Omar |Uma01j 
to estimate the lowest zero of some Dedekind zeta functions. For a fixed support [—X, X] , 
there is a canonical way of choosing a "best" test function, by a method developed for 
the Selberg trace formula in [BSQSj . In order to use the method, which depends crucially 
on a positivity argument, it is necessary to assume the Riemann hypothesis for our given 
L-function. With that caveat, we recall briefly the construction from |BS05j . 

For to G M, let C{X,to) be the class of functions h as above, with the corresponding g 
supported in [— X, X], and the additional restrictions h(t) > for t G M and h(to) = 1. 
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Define 



(3.3) 




Then as X — oo, Fx tends pointwise to the characteristic function of the zeros. Moreover, 
if Fx{to) < 1 for any value of X then to cannot be the ordinate of a zero. Thus, by 
evaluating Fx we can find provable intervals in which the zeros must lie. 

Although the definition of Fx is abstract, it is easy to construct concrete families of 
functions that closely approximate any desired function. For instance, let M be a large 
integer, 6 = X/2M and set 



for arbitrary real numbers a„, 6„. (For self-dual L-functions, we restrict to even test func- 
tions, i.e. all bn = 0, and divide the final formula by 2.) On the other side of the Fourier 
transform, this corresponds to taking g = f * f, where / linearly interpolates arbitrary 
values at multiples of 6. 

The sum over zeros in ()3.3j) is then a positive definite quadratic form in the numbers 
a„ and 6„. To compute the matrix of the form essentially involves computing the explicit 
formula for functions g that are translates of a fixed function of small compact support. 
That requires almost no extra work, since we may compute the formula for all localized 
test functions simultaneously. Once the matrix is known, the infemum in ()3.3p over this 
restricted class of test functions is easily found as the minimum of the quadratic form 
subject to the linear constraint /i(to) = 1- This involves inverting the matrix, after which 
the minimum may be found quickly for many different values of to- 

For an L-function of degree r and conductor A^, the density of zeros at height T is roughly 
^logA^(^) . Therefore, in order to resolve features around height T, the uncertainty 
principle says we should know the numbers c„ for n up to about N(^^y. In the self- 
dual case, the extra division by 2 replaces this by its square root; thus, the complexity 
is on par with that of the approximate functional equation or the algorithm of Section 
EJ although it is much more sensitive to the local spacing of zeros. (Heuristic arguments 
based on experiments and random matrix theory |Odl87j indicate that the minimum gap 
between zeros can be arbitrarily small relative to the mean value; although such small 
gaps are expected to be very rare, we could in principle need many more coefficients than 
for the "typical" zero at height T.) In practice, the explicit formula is clean and easy to 
implement since there are no error terms to estimate with functions of compact support. It 
is particularly well-suited to finding low zeros or to situations where the numbers Cn may 
be computed quickly, as is the case for Artin L-functions; cf. Section 16.11 

As mentioned above, the minimization procedure requires assuming the Riemann hy- 
pothesis. If one is willing to do so, the method may be made completely rigorous, and may 
even be used in place of Turing's method for verifying Artin's conjecture. However, it is 



(3.4) 
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more natural to use it as a quick check in order to fine tune and validate the subsequent 
rigorous methods. In fact, it is helpful to assume Artin's conjecture and apply the method 
to the irreducible Artin L-functions directly. That thins out the spectrum, making it easier 
to isolate individual zeros. We have carried out this procedure for a few examples in Section 



4. Turing's method 

Turing's method for verifying the Riemann hypothesis is described well in his paper 
jTur53j , although there are some errors in the details that were later corrected by Lehman 
|Leh70j . The method has subsequently been extended to Dirichlet L-functions by Rumely 
|Rum93j and Dedekind zeta functions by Tollis |Tol97p . Our contribution is to work out 
the details necessary to apply it to an arbitrary L-function with simple zeros. 

Our argument essentially follows that of Turing. We begin by setting some notation to 
be used only in this section. For t not the ordinate of a zero or pole of A, let 

1 /"^/^ r' 

(4.1) S{t) := / —{a + it)da. 

By convention, we make S{t) upper semi-continuous, i.e. when t is the ordinate of zero or 
pole, we define S(t) = lim£_,o+ S(t + e). 

Next, for ti < t2 let N{ti,t2) denote the net number of zeros with imaginary part in 
(^1,^2], counting multiplicity. When neither ti nor t2 is the ordinate of a zero or pole, we 
may calculate N(ti,t2) using the argument principle, as in ()2.ip . Let C be the rectangle 
with corners at 2 + iti, 2 + it2, —l + it2, —l + iti, with counter-clockwise orientation, H 
the half plane |s G C : 3?(s) > |}, and its complement. Note that by the functional 

equation, we have j^{s) = ~ Hence, 

= —53 / —{s)ds = —'^ i — (s)(isH — $5 / — [s)ds. 
^ JcnH A JcnH 7 ^ JcnH L 

Now for the integral of L'/L we move the right edge of the contour out to cxd, where the 
integrand vanishes. We thus obtain 

(4.3) N{tr,t2) = -531og7(s)|!^;;; + S{t2) - S{t,). 



^Tollis applied his method to cubic and quartic fields. In these cases, there is a slight advantage in 
passing to the normal closure and separating into irreducible Artin L-functions, as we have done for the 
A5 cases in Sectional 
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We select a particular branch of log 7(5) by using the principal branch of logF. With 
this choice, set 



(4.4) $(t):= 



and 



71 



\ogN logTT/ , , c-v;^ \ , c-\^i r 
arg e H — t — | rr + ^ y, /^j I + ^ / , ^^S ■'- 

j=i J j=i 



N{t) := $(t) + 5(t). 



1/2 + it + Hj 



(4.5) 

Then N{ti,t2 
i then N{t) ■ 

standard reference point, so only changes in N{t) are meaningful. (Put another way, the 
branch of log 7 chosen in (|4.4j) is non-canonical.) For large t, $(t) may be evaluated quickly 
by an effective version of Stirling's formula: 



= N{t2) — N{ti). Note that if L is self-dual and vanishes to order < 1 at 
N{0,t). In the general case, although we still have N{t) G Z, there is no 



(4.6) 



3 log r (2;) 



53 



+ e 



8\'^{z) 



for 2 e c \ : 



where the notation / = Q{g) means |/| < g. 

Turing's method is as follows. Recall that A(| + it) is real valued. Thus, if we have 
an accurate procedure to compute A(s) then we may locate all simple zeros on the line 
9ft(s) = I by observing its sign changes. If it turns out that all of the zeros between 
ordinates ti and ^2 are simple and on the line, then we can deduce the Riemann hypothesis 
in that interval by computing N{ti,t2) (minus the contribution from any poles between ti 
and ^2) and finding the same number of sign changes over the interval. 

To compute N{ti,t2), we could evaluate ()4.2|) numerically. However, this would require 
many evaluations of A(s) and would be difficult to carry out rigorously. Fortunately, Turing 
devised a simpler method, based on the fact (first due to Littlewood for ({s)) that S(t) 
has mean value 0. Thus, the graph of N{to,t) — $(t) for any fixed to oscillates around a 
constant value; if we were to plot the same function using the measured number of zeros in 
(tQ,t], then any zeros that we had missed would be obvious as jumps in the graph. 

This can be made rigorous as follows. Let to be a large number that is not the ordinate 
of a zero or pole, and assume that between ordinates to — h and to + h (for some h > 0), 
we have located several zeros of A(s), i.e. we have found small intervals (a„,6„) such that 
A(i-|-mn) and A(|-|-i6„) have opposite sign. Let A^ieft(to,t) (resp. Alright (^0, ^)) be the step 
function which is upper semi-continuous, increases by 1 at each a„ (resp. 6„) and vanishes 
at t = to- We then have 
(4.7) 

N{t) < N{to) + N^tito, t) for t < to and N{t) > N{to) + N^^uito, t) for t > to. 
From these, we can deduce upper and lower bounds for A^(to); integrating ()4.7|) . we get 



(4.8) N{to)h + 



to+h 



Nnght(to,t) dt < 



to 



to+h 



N{t) dt 



to 



to+h 



<l>(t) dt + 



to 



to+h 



Sit) dt 



to 
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and 
(4.9) 



to 



N{to)h+ / NMt{to,t)dt> 



to 



to 



to 



N{t)dt= I ^{t)dt+ I S{t)dt. 

If we have in fact located all zeros in the interval (to — ^, + h) with some amount 
of precision (as measured by the size of the intervals (a„,6„)), then we can expect these 
bounds to be close to the truth. Moreover, if we have effective upper and lower bounds for 
the integral of S{t), then for h large enough, ()4.8|1 and ()4.9|1 will bound a single integer, i.e. 
we can unambiguously determine N(to). Doing this for two different values to = ti,t2, we 
obtain N(ti,t2). 

One nice feature of Turing's method is that precise knowledge of the zeros is only required 
in the short intervals around ti and t2, and even there one can make a trade-off between 
the precision of the zeros and the length h of the interval. For the bulk of the zeros between 
ti and ^2 it suffices to observe the sign changes. 

The remainder of this section is devoted to bounding J S{t)dt, cf. Theorem 14.51 below. 
Our starting point is the following formula, obtained by Littlewood's box principle (see 
fnm §9.9]): 

"ti 



(4.10) 



S{t) dt 



1/2 



log \L{a + it2)\ da — / log \L{a + iti)\ da. 



1/2 



Lemma 4.1. Let notation be as above, and set B := supsj{(s)=3/2 
strip [seC: -\< 3?(s) < §}, 

P{s + lfP{s-2 



L(s)p. Then, for s in the 



(4.11) 



\L{s)\' < B\x{s)Q{s)\ 



PisYPis - 1 



Remark. The power of |(5(s)| in the above is not optimal; for 3?(s) = ^, the "convexity 
bound" says that we can put instead |<5(s)|^/^^^, with a constant depending on e (see 
while the Lindelof hypothesis would have |(5(s)|^. Our present choice permits us 



to avoid Stirling's formula in the proof, and thus obtain a clean bound that is uniform in 
all parameters. 

Proof. We consider first the case when L{s) is entire. Set 
(4.12) F{.s) := L{s)L{l - s) = xisy'Lisf. 

Plugging in the definition of x{s), 

7(0- + it) 



(4.13) 



\F{a + it)\ = \L{a + it)\' 



L{a + it)\' 



7(0- + it) 



7(1 -a + it) 



7(1 — a — it) 

Note that when a = | + a positive integer, the ratio of 7 factors reduces to a polynomial; 
in particular, 

(4.14) \F{3/2 + it)\ = \L{3/2 + it)\^\Q{-l/2 + it)\ < B\Q{3/2 + it)\. 
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The inequality holds since 3?(/ij) > — | for all j. Next, from the functional equation we 
have F{s) = F{1 — s), so that 



(4.15) 

Hence, by (gllll), 
(4.16) 



\F{a + it)\ = |F(1 -a + it)\. 



|F(-l/2 + 2t)| <B\Q{-l/2 + zt)\. 



Thus, the function F{s)/Q{s) is bounded by B on the lines 3fJ(s) = — ^ and 3fJ(s) = |. Note 
that although Q{s) has zeros, F{s) has trivial zeros at the same points; in fact 



(4.17) 



F{s) 

Q{s) 



A(s)L(l 



lis)Qis) 



A(g)L(l - s) 
7(^ + 2) 



Since F has finite order, it follows from the Phragmen-Lindelof theorem that |-F(s)| < 
B\Q{s)\ for all s in the strip. 

If L{s) has poles then the above argument breaks down since F{s)/Q{s) is not holomor- 
phic in the strip. In fact, for each k we get three poles, one at 1 + and two at A^, as 
fl4.17|) shows. To compensate for this, we consider F(s)P(s)^P(s — 1) in the above, in place 
of F{s). One checks that \s^{s - 1)| < |(s + if^s - 2)| on the lines ^{s) = -1/2 and 
3ft(s) = 3/2, so that 

(4.18) \Fis)P{sfPis - 1)1 < B\Qis)PisyPis - 1)| < 5|g(s)P(s + ifPis - 2)|. 

Further, the ratio f^^^l^^^^^^p^g^-^ is holomorphic in the strip, so we may proceed as above. 
The lemma follows. 

Lemma 4.2. Suppose that 

(4.19) {t + Q{iXj)Y > {5/2 + ^{ixj)Y + for some X > 5 and all j = l,...,r. 
Then 

1. Fora e [1/2,5/2], 



□ 



(4.20) 



2V2X 
2. Fora e [-1/2,3/2] 



1 ^4/7r2 + l/4 



X2 



<^^{a + tt)-l\o^ 
7 2 



Q 



+ it 



< 



4r 



(4.21) 



3. For all a, 



(4.22) 



^P' , . , maxiam, 0) 

^-ia + zt)< — 
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(4.23) 



^ a + «t = -log— + - > — ' 



We apply the Stirling-type estimate |Leh70j 



(4.24) 

This yields 
(4.25)^ 

3fJ^(a + zt)-^log 
7 2 



r' 1 

-(z) = logz-- 







2/n'- 



1^(^)2 -gfj(z)2 1 



for 3fJ(z) > 0. 



Q 



+ it 



1 



E 



log 



3/2 + it + fij 



a + it + Hj 



a + it + fij 



+ e 



|(t + 53(/i,))2-(a + 3?(/i,))^ 



For the lower bound in ()4.20|) we need an upper bound for the expression in brackets. 
By hypothesis, the 9 term is bounded by For the others, put a + it + fij = x + iy, 

(5 = 3/2 — a, so that \P\ < 1 and x and y are constrained by y^ > x^ + X^, a; > 0. Then, 
using the inequality log(l + u) < u, we have 



log 



(4.26) 



X + P + iy 



X + iy 



+ 3?- 



X + iy 



.log 1 + 



2/3x + 



x^ 



y 



+ 



X 



x^ + y^ 



^ (/j+l)x + /3V2 ^ 2a; + l/2 ^ 1 ^ 1 



2X2' 



There are r such terms, and the lower bound follows after multiplying by — |. 

The upper bound is similar, but uses the second order inequality log(l + u) >u — /2. 
We omit the details. 

2. Similarly, 

(4.27) 



< 



a + 3?(yUj 



< 



^ ^ + + -^2(a + 3fJ(/.,))2 + X2 - y2X' 



3. 

(4.28) 



.P' 



^—{a + it) = 



k=l 



a + it — \k 



< 



k=l 



max(o", 0) ^ max(crm, 0) 



(t2 + X2 



X2 



(Note that this bound is of faster decay than estimates 1 and 2. That is because we 
have control over the real parts of the poles, while nothing prevents ^{fij) from being of 
comparable size to X. One could also obtain an 0(X~2) bound in 1 and 2, with a constant 
depending on the fij.) □ 
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Lemma 4.3. Let w be a complex number with \^{w)\ < |. Then 

onN f\ {x + l + w){x + l-w) , ^ ,,^f 1 1 \ 

4.29 / log ^ — dx < (log Am —— + = . 

Jq [X + w)[x — w) \1 + w I — w J 

Proof (sketch). Note first that equality is attained at w = 0. Using the principal branch of 
the logarithm, set 

(4.30) f{w) = log ^ ^ dx and g{w) = + 



^0 V {x + w){x — w) J 1 + w 1 — w 

These define analytic functions on C \ M, with real parts extending continuously to M. 
Further, ()4.29j) is equivalent to the assertion 

(4.31) m{w)-i\og4)g{w))<0 for |3i(«;)| < ^. 

Note that f{w) and g{w) are each asymptotic to as \w\ oo. Thus, ()4.3H] holds 

for Q{w) sufficiently large; we check that in fact |Q'(w)| > 2 is enough. By symmetry and the 
maximum modulus principle applied to the function e^^'^^~^^°^^^^^'^'' on the rectangle with 
corners at ±| and ±| + 2z, it suffices to check (j4.31|) on the real axis and for ^{w) = |. On 
the real axis we calculate the integral explicitly and verify the inequality using calculus. For 
5R(w) = I the inequality is strict, so we may verify it computationally for < '^{w) < 2. □ 

Lemma 4.4. For a > 6 + 1, define 
Then 



z; 



.V Z' 



(4.33) r^{a)>^-{a + it)>r^{a) 
and 

(4.34) zeiaY < \L{(j + it)\ < ZeiaY . 
Proof. From (jl.ip we have 

(4.35) r|^(a) - 5?|(a + ^t) = J^p'^^ logp UEK.^'")' " \{v'' ^ V"'') • 

^ p k=\ \ i=i / 

Pairing the terms for and a^^, we see that each summand is < 0, from which the 
second inequality of (|4.33p follows. The first inequality is similar. For (j4.34p . integrate 
flOH|l from a to oo. □ 

Theorem 4.5. Suppose ti and t2 satisfy ()4.19|) . and set 

(4.36) Ce:=\ogZe(Pj+ Tlog^da- log z,(a) rfa + (log4)^ f 

V2/ J 3/2 J 3/2 Ze\2J 
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where ze and Zq are as in Lemma \J^\ (In particular, cq ~ 5.65055 J Then 
(4.37) 



-K j^' S{t)dt<Uog Q{^ + it, 



+ 1ok2 



log 



Q 



+ iU 



+ cer + 



v^(X - 5) ■ 



Remark. Note that there is no assumption on the order of ti and so one obtains a lower 
bound as well by reversing their roles. 

Proof. By ()4.10|1 . we need upper and lower bounds for j'^^ log \ L{a + it) \ da. For the upper 

bound, we use Lemma (4. II 

(4.38) 

r'i/2 POO 

log |L(cr + zt)| (i(T = / log \L{a + it)\ da + / log \L{a + it)\ da 



1/2 



1/2 



< - log 5 + - 
- 2 ^ 2 



1 /"^/^ 



log 



1/2 



J 3/2 

7(1 ~ a + it)Q{a + it) 



7(cr + it) 



da 



1 

+ 2 



3/2 



log 



1/2 



P{a + l + it)^P{a -2 + tt) 



P{a + it)^P{a -1 + it) 



da + I log \L{a + it) \ da. 

3/2 



We bound the first and last terms with Lemma [4.41 



(4.39) 



5R(s)=3/2 



sup |L(s)| + / log \L{a + it)\ da < log Ze 



3/2 



log Zg{a) da. 



3/2 



For the second term, we replace a by 2 — cr in the top 7 factor, and use the recurrence for 

7 to get 



(4.40) 



3/2 



log 



7(a + 1 + it) Q{a + ii) 



da. 



1/2 7(0- + it) Q{a-l + it) 

By the mean value theorem, the integral equals 3? [^(cr* + it) + ^(o"* — 1 + it)'\ for some 
a* e [1/2,5/2]. Thus, by Lemma 1^ dOHll is at most 



(4.41) 



log 



Q{- + it 



vr2X2 2V2X 



Similarly, we see that the third term is bounded by 

We now turn to the lower bound. This part is more delicate since we must take into 
account the contribution of zeros near | + it. We use Turing's idea of comparing log \L{s)\ 
to log |L(s + l)|; the difference between these looks like a value of the logarithmic derivative, 
which we can make precise with the help of Lemma 14.31 

Proceeding, we first clear the poles of A by writing F{s) = A{s)P{s)P{s — 1). This 
function then has the Weierstrass-Hadamard product 



(4.42) 



F{s) 



- I e 



s/p 
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where p runs over the zeros of A and 3?(a) = — X]p^(p)- Next, we spht the integral as 
follows: 
(4.43) 



|.3/2 

\og\L{a + it)\ do = I log 

1/2 Jl/2 

P{a + l + It) 



F{a + it) 



f3/2 

+ / log 

'1/2 



P{a-1+ it) 



F{a + 1 + it) 

5/2 



■3/2 

do + I log 

1/2 



7(cT + 1 + it) 



da 



da + I log \L{a + it) \ da 

3/2 



7(cr + it) 
log \ L{a + it) \ da. 



3/2 



The second term may be estimated, as above, by Lemma f4. 21 and the mean value theorem: 



/•3/2 

(4.44) / log 

Jl/2 



j{a + l+it) 



7(0- + it) 



da > - log 



Q[- + it 



1 ^4/7r2 + l/4 



X2 



The third term is positive since ^ 1 for cr > 0. The fourth and fifth terms are 

handled by Lemma [4.41 

As for the first term, from (j4.42j) we have 



(4.45) 

Thus, 
(4.46) 



log 



F{s) 



F(s + 1) 



X _ £+1 
P 



s + 1 - p 



s- p 



3/2 



log 



1/2 



F{a + it) 



F(a + it + 1) 



da = -y2 I fog 

'1/2 



3/2 



0" + it + 1 — p 



a + it — p 



da. 



Now, by the functional equation, the zeros of A either lie on the line 3fJ(s) = | or come in 
pairs p,l — p. Applying Lemma 14.31 with w = ^ + it — p, we see that ()4.46p is bounded 
below by 



(4.47) 



it — p 



Again by ()4.42|) . this equals 
(4.48) 

-{\og4)^^(- + it 



F \2 



= -(log4)3? 
> -(log 4) 



— - + it] H + it] H + it] H + it 

7V2 J L\2 J P\2 J P\2 

4r 4 /3\ 2m" 
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Altogether, we get 



j^hg\L{a + it)\da> Q-log2^ log q{^ + it 



(4.49) 



/>5/2 /"OO 

/ \ogZ0{a)da + / log2;e(cr) (ia - (log4) 

J3/2 J3/2 

1 ^ j,\og{Ae) + \\ 4mlog2 



4/3 
^A2 



Finally, we combine the upper bound for t = t2 and lower bound for t = ti. We get the 
stated main term plus error 



(4.50) 



V2X 



^log(16e3)+|\ 4(log2) + 5/2 
+ 777^ + m — < 



X2 



y2(X-5)' 



□ 



5. Rigorous Computation of L-functions 

The methods of Section HI depend on a fast, rigorous algorithm for evaluating A(s). We 
describe one such algorithm, based on the Fast Fourier Transform, in this section. We note 
that in the case of the Riemann zeta function, a similar technique was developed and used 
by Odlyzko and Schonhage [OS88J. 

Some algorithms for computing general L-functions were described by Dokchitser |Dokn4j 
and Rubinstein |R,ubn5j . They ultimately boil down to the Cauchy integral formula: 



(5.1) 



- 2^ 



7(^)^(^) 

S - So 



where the contour consists of two vertical lines enclosing sq. Writing L{s) as a Dirichlet 
series and using the functional equation, one is lead to study integrals of the form 



(5.2) 



27ri 



7(_s )n 



ds. 



So 



taken along a vertical line far to the right. Rubinstein, following an idea of Lagarias and 
Odlyzko |LU79j . inserts a factor designed to cancel the decay of the 7 factor, e.g. e^~^^ for 
some 7] close to ibl. Without this factor, very high precision is required to calculate A(s) 
when Q^(s) is large. 

These algorithms are good when one is interested in computing A(s) at specific points, 
e.g. for locating zeros of L{s) precisely. They suffer from the disadvantage of being difficult 
to carry out rigorously, basically because ()5.2j) is a two parameter family (indexed by sq 
and n) of integrals, for which uniform asymptotics are hard to obtain in certain transition 
ranges. 
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For Turing's method, we need an algorithm for rigorously computing A(s) for many values 
of s, not necessarily with high precision. For that we consider instead the one parameter 
integrals 

(5.3) [ A{s)e-'' ds and [ -f{s)e-'' ds. 

2m J zm J 

These are essentially Fourier transforms, and they contain enough information for evaluat- 
ing A(s) quickly, if one is interested in many points. They also involve only a single Mellin 
transform, making rigorous computation more accessible. 

Precisely, let ri G (—1,1) and set F(t) := A(^^ + itj e~^^ . Then the (inverse) Fourier 
transform of F is 
(5.4) 

F(x) := r F(t)e~'"* dt = ^ f A(s)e("+*'?'')(^/2-^) ds 

271" J-oo 

= — / A(s)e(^+^^'')(i/^-^) - V Res.=,A(s)e(^+^^'7)(V2-.)_ 

The residue sum is straightforward to evaluate assuming we have complete information on 
any poles of L{s). We multiply the Euler product ()1.1|) out to a Dirichlet series, writing 
-^(■5) = Yl'^=i ^n^""*- Then the first term of (j5.4j) is 



00 „ 00 ^ 

(5.5) Ya^— / ^(s)e(^+^^^)(^/2-V^rf, = ey ^G(x + log^ 



where 



(5.6) Giu; V, {^,}) := ^ / e^^+^^^^^^Z^-^) fl T^{s + /i,) ds. 

Let us assume for now that we have a procedure to compute G{u] rj, {yUj}), and thereby 
F[x), to prescribed precision; we return to this point in Section [5.11 below. In order to 
use the FFT to compute F from F, we first need to discretize the problem. To that end, 
let A, B > he parameters such that q = AB is an integer. By the Poisson summation 
formula, 
(5.7) 

E^/^ , ^\ 271 f 2Tik\ f km\ 27i \ ^ f mn\ sr-^ f 2'Kn 

^(i + *^) = iFE^(ir)<iij=iF E <t)^H^ 

fceZ fceZ ^ ^ ^ ^ n(modg) ^ ^ ^ fceZ ^ ^ 

Thus, the functions F(m) := Efeez^(f + and F{n) := Efcez^(^ + 27rAA;), which 
are periodic in m, n with period g, form a discrete Fourier transform pair. 
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Note that since F is real-valued, F{—x) = F{x). Thus, for \n\ < q/2 we have 

'27m 



(5.8) F{n) = F 



^7m 



B 



oo 
k=l 



B 



^ k=\ ^ 



For A even moderately large, the terms for > 1 fall within the asymptotic range. Precise 
bounds are given in Section [5.31 below: in particular, we may apply Lemma 15.41 with x = 
27i{A ± ^) to compute the sums over k. Hence, it suffices to calculate F{2Trn/B) for 
< n < q/2. On the other hand, to compute F{m/A), we need to bound the terms of 
F{m) for 7^ 0. We have already obtained a suitable bound for the L-function in Lemma 
14.11 The sum of this bound over 7^ is the content of Lemma 15.51 

5.1. Computing 77, {/Zj}). For brevity, some of the results of this section are only 
sketched. Our emphasis is on the details necessary for rigorous computation. For more 
general background information we refer the reader to [Boo03a| IDok04| IRub05j . 

One simple method for calculating integrals such as (j5.6|) that is easy to make rigorous 
is the power and log series, obtained by shifting the contour of (|5.6|) to the left: 



(5.9) 



G{u] r], = ^ Res, 

pec 



poles p 



where P{u; p, 77, {pj}) is a polynomial of degree one less than the order of the pole at p. 

For example, in the case of Galois representations, the pj are all either or 1, and the 
residues in ()5.9j) may be evaluated by the following: 
(5.10) 



-2n) 



s + 2k {2k)\\ 



exp 



^2n 2 



g(7re^) (s + 2k) 



00 



+5:i((-i)^2-cw+E|^ 



-27t) 



(2A;- 1)!! 



exp 



E 

V71=l 

OD 



i=2 

1 




log(47re'^ 

2n-l 2 




1) 



j=2 ^ \ n=l 



^ (2n - ly 

n=l ^ ' 



(S + 2A;- 1) 



for any integer > 0. For general /ij, we need an algorithm to calculate the values of F 
and its derivatives at an arbitrary point in the complex plane; we assume without further 
comment that this is available when necessary. 
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We say that pj and /i^ are equivalent if pj — /i^ G 2Z. For /i ranging over an equivalence 
class, the functions T^{s + fi) share all but finitely many poles. Thus, (j5.9p may be broken 
naturally into parts corresponding to each class. We can bound the tail of each part as 
follows. 



Lemma 5.1. Let p be a pole of g{s) 



^(u+i^rj){l/2-s) 



YYj^iTR^s + pj) of order n, with 



Let Cj be 



3fJ(p + P,) < for J = l,...,r and (27r)'^e2« < i n;=i(|2 - p - /i, | - 1) 
the coefficients of the polar part of g around p, i.e. such that g{s + p) — Yl^=i'^j^~ 
holomorphic at s = 0. Then 

oo 

(5.11) y^ReSs=p^2fc^(-s) 
Proof. First note that 

g{s + p~2) = 



IS 



k=l 



< max Co 



(27r 



(5.12) 



fe'u+^^^g^s + p)l[{s + p-2 + p,: 

(-27r)V2«+^¥^ g{s + p) 



-1 



n;=,(2-p-ft) n;-.,(i 



:) 



Next, let f{s) be a meromorphic function with polar part ais~^ + . . . + a„s~" at 0. If x is 
a complex number with |x| < 1 then the function ^^f^ has polar part a'^s~^ + . . . + a'^s"", 
where a'- 



ELo ^i+fc^;''- Thus, 



(5.13) 



max \ a!j\ < 



max 




1 - 


\x\ 



Let c']^s + . . . + c^s 
see from (|5.12p that 



max I c'- 1 < 



be the polar part of g{s + p — 2). Applying (j5.13|) r times, we 
(27r)^e2" 



max I Cj I 



(5.14) 



n 



(27r 



|2-p-/ij| 



) 



\r ^2u 



max Co 



n 



p-Pj\ 



< 



max c 



-J I 



Repeating this procedure, we see that the coefficients of the polar part of g{s + p — 2k) are 
< 2"^^ max |cj|. The conclusion follows. □ 

The lemma says roughly that if we compute the residue sum for all poles with real part 
down to 3?(p), the tail of the series (from poles at p — 2k) may be bounded by the data 
from the last term added. Moreover, ()5.12|) gives an algorithm for computing the data at 
p — 2 from that at p, and shows that the terms eventually decrease factorially. Thus, we 
may use this to compute G{u] rj, {pj}) to any desired precision for a given u. 
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Since the /ij are arbitrary, this procedure is general enough to compute derivatives of 
G{u] 1], {/Xj}) as welL For instance, for any k we have 



(5.15) 



G\u; f], {fij}) = {jik + 1/2)G(m; ?7, {nj}) - 2nG{u; t], {/x' }). 



where fi'j = fij if j ^ k and ^'j^ = fik + 2. Higher derivatives may be computed in a similar 
fashion. (In fact G satisfies an rth order differential equation, due to the recurrence for F; 
thus, the derivatives of all orders are determined from the first r.) 

Note that for u large this method requires high precision due to cancellation, and is 
therefore inefficient. The essential point that makes it worthwhile is that for a given 7 factor 
the calculations need only be performed once, as one can develop local approximations to 
G{u; 1], {fij}) for later rapid evaluation. The computation may then be recycled and used 
for any L-function with the same fij] this is useful for functions in an arithmetic family, 
such as Artin L-functions. 

More precisely, suppose we wish to calculate G{u] r], {fij}) for u in an interval /. Choose 
e > and sample points Um such that each m G / is contained in a unique interval 
[um — £, Urn + £)■ For u in the mth interval, we have by Taylor's theorem 
(5.16) 

G{u;r],{^j}) = > ^ — . : ' (u - Ur 



K-l 
A:=0 



G^''\u^;r],{ijj}) 



k\ 



+ O I max 

\u*-UTn\<e 



K\ 



We may evaluate the derivatives precisely using ()5.9|) . As for the i^th derivative, a uniform 
bound is obtained by shifting the contour of (j5.6p to 3?(s) = |: 



(5.17) 



\G(^\u■,v,{^^,})\ 
Kl 



1 

< — 



K\ 



n 



^R[^ + 'it + fij 



dt. 



For large K this is of size f ^^(i-??) ) ' ^^^^ iong as e is small compared to ^(1 — r]), 
we may compute and store the coefficients of ()5.1(i|l . yielding a fast method to calculate 
G{u; T], {fij}) for any u E I. 

Moreover, we can improve the efficiency of our algorithm if the sample points of (j5.8|) 
coincide with multiplies of 2e, i.e. if ^ G Z. For any given sample point x, we approximate 
F{x) via a truncated series (keeping track of the error terms from ()5.16|1 and Lemma (5.31 
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below) : 



IVl / 

^ 'ItPV ^ ^' ^^^^ 

n=l ^ ^ 



AT 



(5.18) «E E ^E ^"''^:r"''"'" C°g 



m log-^e[u„-e,«m+e) fc=0 



fc=0 m 



where 



(5-19) Si^^:= Y: ^(^^ 

log ^^e[Um-£,Mm+e) 



k 

n 



Since x + Mm is another sample point Um', the fcth term of ()5.18p is a convolution of the 

sequences (indexed by m) ^"'^"-j"'^^^" and S^r^^; thus, we may evaluate it efficiently for all 
X simultaneously by appeahng again to the FFT. 



5.2. Complexity. We may now consider the complexity of the algorithm. Note that by 
Stirling's formula, F(t) decays roughly like e"*^^"''-'^* for t > 0. Ideally we should choose 
1 — ?7 of size in order to compute values up to height T. Adjusting the constant of 
proportionality (i.e. choosing rj relatively close to or far from 1) allows us to trade off the 
computational precision and number of coefficients needed to overcome the error terms 
below. Finding a good compromise between these two is best done by trial and error; cf. 
Section ini 

Since S x T^^, Lemma f5.HI shows that in order to compute ()5.18|1 we need on the order 
of viVT^ terms of ()5.5p . or roughly the square root of the analytic conductor. Note that 
the values of A and B enter only in the Fourier transforms, and do not significantly affect 
the computation of (|5.18p . We set B equal to a multiple of T, depending on the chosen 
value of f]. As for A, as mentioned in Section |2l the density of zeros of F{t) around height 
T is 2^ log A^(^) ; one can expect to take A equal to a multiple of this. Thus, this method 
has complexity consistent with computing a single value by the approximate functional 
equation, after which we get many values in mean time Oe((A^T)^), which is essentially 
best possible. The gain comes from the fact, as emphasized above, that only a single 
G-function is involved. 



5.3. Asymptotics. To complete our understanding of G{u] rj, {yUj}), in order to accurately 
calculate ()5.5|1 . we need an asymptotic bound for large u. If we write fx = + ^(l + 
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Y7j=i f^j) then by the method of stationary phase, we have 



(5.20) G{u;r],{fi,}) = ^±_e'^("+^^'') exp(-ei("+*^'')) (l + 0(e-2"/'^)) , 

where the imphed constant depends on the iJ,j. For r = 1 the formula is exact, i.e. the O 
term is 0. For r > 1, one can work out exphcit constants case by case, which is preferable 
if sharp error terms are desired. Otherwise, we get a bound that is close to ()5.20j) simply 
by shifting the contour of ()5.6|) to the right. 

Lemma 5.2. Let 5 = f (1 - \v\), = + Yr' = + + E-=if^)' ^ = 



2^2:±i£!l^e-^^, and X = vrr^e-^e^^A. Then for X>r, 
(5.21) \G{u■,r^,{^,,})\<Ke''(^>e-^f[[l + '-^)' 



Remark. This is within a factor 0(5^^/^) of the correct asymptotic if 5 <C re""/''. 



dt. 



Proof. We write s = 2a + 2it in ()5.6|) to get 
(5.22) 

\G{u- V, {/z,})| < vr-'('^+^+i)-^e-(^2'^)-"^ H f[\r(a + zt + ^) e"?(2*+^(^^)) 
Applying Holder's inequality, we get integrals of the form 

/oo 
|r(a + it)fe^*rft, 
■oo 

where a = a + 3?(/ij)/2. Assuming a > 1, we may apply the inequality 
(5.24) |r(a + it)| < v^(a+ |t|)"-^/V"l*l/2 for a > ^. 

(To see this, note that |r(i+2t) | = V Trsech Trt, use the recurrence for F and ()4.24|) to reduce 
to the region |<a<|,0<)f:<2, where the inequality may be checked computationally.) 
Thus, we have 

|F(a + it) Te^* dt < 2(27r)''/2 / (« + tyi'^-^^) q-'^' dt 



, xr./9 Krn^{r{a- 1/2) + 1 

(5.25) < 2(27r)'/2e^- ^ ^ < 



< 2(27r) 



(^^^)r(a-l/2)+l 

+i p^/e^a-l/2 + l/2rV^""'^'^^'^' 



(5r 
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Substituting this bound into ()5.22|1 and collecting terms we obtain 
(5.26) ire«(-)- n (^) ^^'^ < Ke^^^^e- \{ (^) (l + ^ 

The result follows upon taking a = X/r. □ 

With this bound in hand, we can estimate the error in truncating the series ()5.5p . 

Lemma 5.3. Let M be a positive integer, x G M. Let 6, 1'jjfi, K he as in Lemma \5.^ and 

set X = T:r5e~^ {e!^ / \/nY^^ ■ C, a > be such that |a„| < Cn°' for all n, and put 
c = gfj(/i) + l + a, c' = max(cr/2 -1,0). Then for XM'^/'' > max(c', r), 



(5.27) 



n>M 



< 



rll l + TxA 



Remark. Different values of C and a are appropriate for different ranges. For small M, one 
can take C = 1, a = logg r + 6, while for larger M it is better to choose a smaller value of 
a and compute C from the coefficients. 

Proof. Using Lemma [5.21 we have 



(5.28) <k'y: 

n>M V \ ViV / 



where 



The condition on X ensures that the terms of ()5.28j) are monotonically decreasing. Thus, 
we can estimate by the integral 



(5.30) K' / f "^e-^*'"- rft = fll(X'^-/Y / y'^'-Z^-^e-^ < fll^^ . 

^ ' Jm 2 ^ ^ JxMVr 2 XM2/'--c' 



□ 



The next two lemmas bound the error introduced in discretization. 
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Lemma 5.4. Let x G M, A > and let notation be as in Lemma 1X31 Then for X > 



maxl^c ,r), 



'^F{x + 27ikA) 



(5.31) 



k=0 



Ress=p _ ^2-kA{i/2-s) 

pe{l+Afc:l<fe<m} 



+ e 



K 



1 - e- 



-X 



1 + 



Cr/2 
X-d 



Proof. The residue sum comes from summing the polar part of ()5.4|) with x + in place 
of X. For the rest, we apply Lemmas 15.21 and 15.31 (with M = 1) to get the bound 



(5.32) 



K 



-X 



Cr/2 



X-d 



rvj \ 



X 



for the A; = term. To pass from this to the /cth term, we multiply by a factor not exceeding 



(5.33) 



g2.fcAJ?(^)g-X(cxp(4.fcA/r)-l) ^ g^p ( 27rfcA 



2X X 



< exp -2'nkA 



r 2nkA 
1 2 4:7ikA 



ATTkA\ 



«+---+ 

2 r r 



< e 



-TrfeA 



The result follows on summing the geometric series. 
Lemma 5.5. Let t eM. and put s = ^ + it, 



□ 



(5.34) 

and 

(5.35) 



E = Zei3/2r\^is)\e'i 



Qis. 



P(s + l)2p(s-2) 



^ vrr 1 X - 

p = — y, arctan 



P(S)2P(S- 1) 

3?(s + /ij) 4 



1/2 



4 2 



E 



3(s + /ij)| 7r2 |5>(s + /ij)2-SR(s + /ij-; 



;i) // $5(s + fij) > for all J = 1, . . . ,r and p - fr] > then 



(5.36) 



J2F{t + kB) 



k=0 



< 



E 



1 - e 



(2) // ^^(s + fXj) < for all j = 1, . . . ,r and p + fr] > then 



(5.37) 



< 



k=0 



E 

1 _ e-(/3+¥^)^' 
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Proof. We treat only the first case, the second being similar. Lemmas 14.11 and 14.41 imply 
the bound < E. We consider the same bound with t replaced by t + kB. Note that 

if \Q{s + fij)\ increases for all j then the factor involving P is non-increasing. For the 7 
and Q factors, by the mean value theorem we have 

Q{s + ikB) ^'^^ 



(5.38) 



log 



7(5 + ikB) 



7(s) 



Q{s) 



7 



2 Q 



for some s* on the line between s and s + ikB. Using ()4.24j) . this is 
(5.39) 

.7 = 1 ^ ^' 



< -kB\ 



nr 1 3?(s* + 

y arctan 



$5(s*+/ij)2-gfJ(s*+/ij)2| 
4 



Thus, \F{t + kB)\< Ee-^f^-^"^^^^ . The conclusion follows. 



□ 



6. Numerical results 

We have applied the methods described in the previous sections to a few examples of 
splitting fields of polynomials with Galois group 5*5 and A^, as hsted in Table lUTTl For the 
A^ cases, the Artin conjecture is true for all representations by known cases of functoriality 
|Kim94irTMnT] . That speeds up the process, since we may apply Turing's method to the 
Artin L-functions directly. For the 5*5 examples we verify both conjectures. As expected, 
we found no counterexamples to either conjecture in the tested range \t\ < 100. 



polynomial 


group 


splitting field discriminant 


— 68x — 68 




296360-1^796 


— x'^ — 8x^ + lOx^ — X — 5 


^5 


2 160 396 796 


x^ -x^ + 3x^ - Ux^ - 8x - 8 


^5 


2220 1396 


x^ + 2x^ - 4x^ - 2a: + 4 


A, 


2907330 


x'' + 20a; + 16 


A, 


29O578 


x^ — x'^ + 8x^ — 6x^ + 14a; — 6 


A5 


290^9330 


x^ - 7x^ - 17x^ + 18x + 73 


A, 


24048730 


x^ + 8x^ + 7x^ + 172x + 53 


A5 


2083^° 



Table 6.1. Tested polynomials 



To illustrate the methods, we discuss in detail the 5*5 field of discriminant 2^^3^°17^^ 
given by the polynomial f{x) = x^ — 68{x + 1). Recall that 5*5 has seven irreducible 
representations. We label them 1, x, p^, = Pi ® Xi Ps; P5 = Ps ® X and pg; where 
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X is the sign character and the subscripts indicate the dimensions. As it will turn out, 
the limiting factor in our computations is the conductor of pg? which in our example is 
36081072 = 2^3'^17^. This is the smallest among the table of 5*5 polynomials given in |KMj : 
since that table is ordered by the conductor of p^, it is likely that smaller examples exist. 
(We note, however, that if one is interested only in verifying some instances of Artin's 
conjecture and not the Riemann hypothesis, the holomorphy of L(s, pg) may be checked 
much more easily; there the limiting factor is the conductor of p4, of which |KMj yields 
examples as small as 1609. We have not pursued this possibility.) 

Note that 1, x ^i^nd pg are monomial representations, so Artin's conjecture is true for those. 
Equation (j6.ip below shows that L(s, p^) and L(s, ps) are holomorphic except possibly at the 
zeros of C{^)- Twisting by X; we see similarly that L{s,p'^) and L{s,p'^) are holomorphic 
away from the zeros of L{s,x)- Moreover, we learn from GAP that the representations 
Pa® Pe, P4 © Pe, Ps © P4 © X p'5 © P4 © 1 are all monomial. Thus, in order to verify the 
holomorphy of L{s, p^j, L{s, P4), L(s, P5) and L(s, pg) it is enough to check that L(s, pe) and 
L{s, P4 © x) are non-vanishing at zeros of ({s) and, similarly, that L(s, pe) and L(s, p4 © 1) 
do not vanish at the zeros of L{s, x)- Applying Turing's method to these functions as well 
as L{s,p5 © 1) and L{s,p'^ © x), we can deduce both Artin's conjecture and the Riemann 
hypothesis (up to the tested height) for all representations. 

In what follows we describe the numerical procedure in detail for L{s, po). First we must 
choose a value of r] to use for the computation of G{u] t], {pj})- One can aim to limit either 
the number of Dirichlet coefficients a„ or the precision required in the computation. Since 
the coefficients are relatively easy to compute in our case (we have 2^^ of them), we try 
for the latter. The largest error comes from Lemma (5.31 with x = 0, and is of size roughly 
M^exp(-7rr5e-^(M/ViV)2A), where ^ = |(1 - |r/|), r = 6, = 36081072 and M = 2^^. 
Examining the local factors at small primes we determine that |a„| < 1.26?t,'°S4243 6^ yielding 
c= l + lo 

§4243 6- (The pj in this case are 0, 0, 0, 1, 1, 1.) This error term should be compared 
to the size of the function being evaluated, which is roughly e"'^''*/^. From Theorem 14.51 we 
find that to apply Turing's method up to height t = 100 we need to be able to compute 
the L-function up to about t = 115. Trying a few values of rj, we find that with rj = 0.98 
the error terms are of size 10~^^, compared to 10~^ for the size of the function. Thus, 
with this choice we should use a precision of at least 14 digits; in fact we carry out most 
computations to 30 digits. 

All computations were performed on a 3GHz PC running Linux. They were divided into 
several steps: 



(1) Computing the Dirichlet coefficients a„; 

(2) Estimating zeros by the explicit formula; 

(3) Computing G^^\um: rj, {pj})\ 

(4) Computing S^^^- 

(5) Computing ^fJ^ and L{s,pq) by FFT; 

(6) Turing's method. 
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To ensure correct results, we used the arbitrary precision interval arithmetic package MPFI 
[RR05j for steps 3 through 7. We discuss the steps in more detail below. 



6.1. Computing a^. First we consider methods of computing the coefficients of the L- 
functions L{s,p) for all irreducible representations p. One way is to express them as ratios 
of products of Hecke L-functions, as given by Brauer's theorem; in fact for 5*5 we may 
express each in terms of Dedekind zeta functions of intermediate fields. Precisely, let k be 
the quadratic extension of Q associated to F = Q(xi) C M = Q(xi,X2) where Xi and 
X2 are distinct roots of /, and E = Q{y) C K where y is a root of the sextic resolvent (a 
formula for which is given in |Dum91j ): then we have 



(6.1) 



In turn, we may compute each of the Dedekind zeta functions using the ideal factorization 
functions built in to PARI |The04j . This facilitates the computation of local factors at 



primes dividing the discriminant, allowing us to avoid a detailed study of the possible 
types of ramification. However, it is not well-suited to working out many coefficients. 

Fortunately, there is a faster method that works for all but finitely many primes. Table 
16. 21 shows the unramified local factors for each representation and conjugacy class (labelled 
by the order of elements in the class), where we write x for p~^. For 5*5 it turns out that 
the Frobenius conjugacy class at p is determined by the number of linear and quadratic 
factors of the reduction / of / modulo p, which may be computed from the degrees of 
gcd(a:^" — x, /(x)) for n = 1,2. That computation requires O(logp) multiphcations and 
additions mod p. Thus, by the prime number theorem, for each L[s,p) we may determine 
the Dirichlet coefficients a„ for n < X in time 0{X) (assuming mod p multiplications and 
additions take bounded time, which is appropriate for numbers of the size that we consider). 
Up to the implied constant, that is best possible. Moreover, the technique is very fast in 
practice; we found that it takes approximately seven hours to compute the local factors for 
all p < 232. 



6.2. Estimating zeros. With our computed coefficients, we readily obtain estimates for 
the low zeros by the method of Section |21 Figure lUTD shows graphs of Fx{t), with X = 
log(232)^ for each irreducible L-function. The spikes correspond to zeros, from which we 
get the estimates for the ordinate of the lowest zero of each function shown in Table 16.31 
note that for ( the estimate agrees with the known value 14.1347251417. . . to within the 
precision of the computation. The increase in density of zeros with the conductor and 
degree is apparent in the graphs. Moreover, as the explicit formula is very sensitive to 
errors in the coefficients, the fact that we see spikes of height 1 for the low zeros indicates 
that our coefficients were computed correctly. Each graph took a few minutes to generate. 
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(1 -X)2(l 
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(1 -x)(l -x^) 
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-x^) 
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P5 
P5 
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-X2)2 
-X2)2 


(1-X)(1-X2)2 
(1-X)(1-X2)2 
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(l+x + x2)(l-x^) 
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1 — X 




1 — X 
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1+x 


1 — X 




1 + x 




pi 


1-x^ 


1 + X + x^ 


+ X + X 


(l + x)(l -x^) 




p\ 

P5 
P5 
P6 


1-x^ 

{l-x){l-x^) 
(1 + x)(l -x^) 

+ 


1 + X + x^ 
1 — x^ 
1 - x^ 
(l-x)(l- 


+ X + X 
- X^) 


(1 -x)(l + x3) 
(1 + X + x2)(l + x^) 
(l-x + x^)(l-x^) 
1 - x*5 





Table 6.2. 5*5 unramified local factors 



p 


conductor 


lowest zero 


P 


conductor 


lowest zero 


1 


1 


14.134725142 


X 


3 


8.039737156 


P4 


4009008 


1.108937765 


Pi 


36081072 


0.5717508665 


P5 


36081072 


1.062064850 


P'5 


12027024 


0.8132800720 


P6 


36081072 


1.376872200 









Table 6.3. Conductor and ordinate of the lowest zero of each irreducible L-function 



6.3. Computing G^^^ ?7, {Pj})- Next we compute local approximations of G{u] r], {yU-,}) 

for u in the interval log log . We evaluate 2^^ Taylor series of 16 terms using (|5.9p : 

with these choices, the error term in ()5.16j) is less than 10^^^. This calculation is the most 
delicate, due to high precision and catastrophic cancellation. Nevertheless, the computa- 
tion time for this stage was only a few hours. The graph of \G{u;ri, {fij})\ is shown in 
Figure | 



6.4. Computing Sm\ ()5.18p and L{s,pq). Now we come to the main part of the compu- 
tation, ()5.18p . Most of the time, approximately twelve hours, was spent computing Sm'' ■ 
Note that if we had not adjusted rj to reduce the precision, this calculation could have 
taken substantially longer. 

Once we have Sm\ the computation of ()5.18p and L{s,pq) is very fast. We choose 
B = |Qg(-'232) ~ 1160. Since this is much larger than t = 115, the errors terms from Lemma 
15.51 are negligible. We choose A = ^ 900, which is about 160 times the expected 
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Figure 6.1. Fx{t) for each irreducible representation 



14 



density ^logA^(^)'^ of zeros around t = 115. Thus, the main Fourier transform is of 2^° 
points, which takes only a few minutes to compute. 

Figure 16.31 shows the graph of Z{t) := A(i + it) /\'y{^ + |) which is the analogue of 
Riemann's Z function. We have superimposed the graph of 25Fx{t) over the same range; 
note the good agreement in location of zeros between the two, which gives evidence that 
our computations are correct. Figure shows Z{t) over the higher range t G [90, 100]. 

6.5. Turing's method. Finally, we apply Turing's method to the computed L-functions. 
At the same time, we verify the "working hypothesis" that we can isolate the zeros of the 
irreducible L-functions. This verification takes only a few seconds. 
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Figure 6.2. \G{u; 0.98, {0, 0, 0, 1, 1, 1}) | 




4 6 8 10 12 

Figure 6.3. Z{t) and 25Fx{t) for small t 
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